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Abstract. Transport properties of a gated nanostructure depend crucially on the 
coupling of its states to the states of electrodes. In the case of a single quantum 
dot the coupling, for a given quantum state, is constant or can be slightly modified 
by additional gating. In this paper we consider a concentric dot-ring nanostructure 
(DRN) and show that its transport properties can be drastically modified due to the 
unique geometry. We calculate the dc current through a DRN in the Coulomb blockade 
regime and show that it can efficiently work as a single electron transistor or a current 
rectifier. In both cases the transport characteristics strongly depends on the details 
of the confinement potential. The calculations are carried out for low and high bias 
regime, the latter being especially interesting in the context of current rectification due 
to fast relaxation processes. 


PACS numbers: 73.23.Hk, 73.21.La, 73.22.-f 



Charge transport 

1. Introduction 


2 


In order to meet growing demand for small scale, low-power consuming devices one has to 
downscale transistors and logic circuits and to work with a small number of carriers. The 
natural limit for lowering carrier density is single charge electronics, where phenomena 
such as electric current can be controlled with single electron precision [1,2]. Contrary 
to modern mass production electronics, single and a few electron devices exhibit purely 
quantum mechanical effects such as resonant tunneling [1-5] or quantum entanglement 
[6-8]. Apart from direct applications in nanoelectronics they are also perfect tools 
for probing fundamental problems in single and many-body physics. During the last 
decade a lot of research has been devoted to study the electronic properties of qauntum 
dots (QD) [5,9-14], For sufficiently low temperatures the discretness of the energy 
spectrum of these systems can be clearly visible in transport experiments as Coulomb 
peaks [1,4,15,16] that demonstrate succesive charging and discharging of the QD by 
single electrons. The Coulomb blockade phenomenon, where the charging energy forbids 
an electron to jump to the QD, is the basis for most of the applications of QDs [17-19]. 
Quantum dots arranged into double, triple or more complex systems [20-28] exhibit 
abundance of quantum states which manifest themselves, e.g., in Pauli spin blockade 
[29,30] current or heat rectification [13] effects. 

Another class of interesting quantum systems are quantum rings (QR) [31-35]. Due 
to different from QDs geometry the phenomena observed in QRs are very sensitive to 
phase coherence of the electronic wave function. These are, e.g., the Aharonov-Bohm 
effect demonstrating modification of the electron wave function by a vector potential 
[36] or persistent currents, i.e., ground state currents that flow in QR even without an 
external magnetic field [37-39]. 

In this paper we focus on a complex nanostructure that combines the two mentioned 
above, topologically different components: a quantum dot and a quantum ring. The 
constituents are aligned concentrically (QD is surrounded by QR) so that the system 
conserves the circular symmetry. The dot-ring nanostructure (DRN) has already been 
fabricated by pulsed droplet epitaxy [40,41] with full control of the growth process. It 
can also be made by using atomic force microscope to locally oxidize the surface of a 
sample [12] or by lithography. Another method would be to grow a core-shell nanowire 
[42,43] of, e.g., (In,Ga)As, where the core and shell parts are separated by a tunneling 
barrier. Then by cutting a slice of it one can form a DRN. A DRN exhibits a large 
variety of quantum states [44] what leads to many interesting features. 

It has been recently shown [45,46] that many measurable properties of a DRN, 
like spin relaxation or optical absorption, can be widely changed by a modification of 
the confinement potential of the DRN demonstrating its very high controllability and 
flexibility. These characteristics are mostly determined by the relative distribution of 
the wave functions in a DRN that, in turn, can be changed by external gates or fields. 
The purpose of this study is to demonstrate that also conducting properties of a DRN 
are very sensitive to the details of confinement and that we can realize different single 



Charge transport 


3 


electron devices on this complex structure. Unlike field effect transistors, single electron 
devices are based on intrinsically quantum phenomenon, namely the tunnel effect. In 
this case transport properties are mostly determined by the tunneling rates Us, which 
depend on the overlap of the DRN states with the states of the electrodes. These 
parameters, in turn, depend crucially on the localization of the electron wave function: 
states localized in QD (QR) are weakly (strongly) coupled to the electrodes. Thus Us 
may strongly depend on the quantum state and have to be determined for each state 
individually. This property demonstrates one of the advantages of the DRN over QDs, 
where the possible changes of the couplings are orders of magnitude smaller. 

We discuss charge transport through a DRN in the Coulomb blockade regime near 
the N — 0 -H- 1 transition, i.e., when only one electron at a time can tunnel through 
a DRN between the source (S) and drain (D) electrodes. Throughout this paper we 
assume the magnetic field B = 0 and therefore neglect the electron spin. We show that 
one can tune the device parameters so that it can work as: ( i) single electron transistor 
[4,47] and (ii) electrical current rectifier. The paper is organized as follows: In Sec. 2 
we present a general theoretical background that will be needed to study the transport 
properties of DRNs. In Sec. 3 we demonstrate how by changing the parameters of 
the confinement potential we can control single electron tunneling and build a single 
electron transistor (SET). In Sec. 4 we demonstrate that a DRN can be used as a 
current rectifier. The results are summarized in Sec. 5. 

2. Basic formulas and mechanisms 

We consider a quasi-two-dimensional circularly symmetric dot-ring nanostructure. The 
DRN, placed in the xy plane, is defined by a specific confinement potential, that we 
will discuss below. We assume that the confinement in the z (growth) direction is much 
stronger than the lateral confinement and consequently, the xy-plane motion and the 
vertical one can be decoupled. Then, we can write the electron wave function as a 
product 

il>{r) = i!)\\(r)i/> z (z), (1) 

where the vector r lies in the xy-plane. Additionally, we assume that the electron is 
always in the lowest energy state of a quantum well in z direction and that the potential 
of the well in z direction is infinite. With these assumptions i}) z {z) is given by 

Uz) = /fcos(^) . (2) 

where d is the height of the structure. 

In order to discuss the in-plane confinement potential that forms the DRN we 
introduce the following notation |r| = |(x, y)\ = r. Then the DRN is defined by a 
potential V(r) and occupied by a single electron. The DRN is composed of a QD 
surrounded by a QR and separated from the ring by a potential barrier Vo (r). A cross 
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section and a top view of a DRN with explanation of symbols used throughout the text 
is presented in figure 1. 



Figure 1 . The cross section and the top view of the potential forming the DRN with 
marked bottom of the QD potential (Vqd), bottom of the QR potential (Vqr), top of 
the barrier potential (Vo) and the value of the potential outside the DRN (Vi). We 
assume the bottom of the quantum ring part as a reference value (Vqr = 0). This 
figure shows a DRN with the bottom of the QD below the bottom of the QR (Vqd < 0), 
but also the opposite situation (Vqd > 0) is possible. 


In particular, we assume in our model calculations the radius of the DRN r 0 = 70 
nm. The depth of the quantum well forming the DRN is V± = 90 meV and the zero 
potential energy is set at the level of Vqr, i.e., the potential well offset is equal Vqd- 
The calculations are performed for InGaAs systems (with the effective electron mass 
m* = 0.067 m e ). The results are presented for Vo = 20 meV and for the sample thickness 
d — 5 nm, if not stated otherwise. 

We solve numerically the Schrodinger equation assuming the Gaussian form of 
V(r) [45]. The energy spectrum consists of a set of discrete states E nl due to radial 
motion with radial quantum numbers n = 0,1, 2,..., and rotational motion with angular 
momentum quantum numbers / = 0, ±1, ±2 .... The energy spectrum as a function of 
Vqd is shown in figure 2. The states situated in QD exhibit an increase of the energy 
with increasing Vqd, whereas those situated in QR have the energy (nearly) constant. 
The single particle wave function in the xy- plane, i.e., the DRN plane, is of the form 


V’nO) = T ni (r,0) = R nl (r)exp (il<f >), 


(3) 
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Vqd [meV] 


Figure 2. Energy spectrum of a DRN as a function of the position of the bottom 
of the QD potential Vqd- The circles indicate states for which wave functions are 
presented in figure 3. The dotted line shows the energy of the lowest two-electron 
state (see text in Sec. 3). 


with the radial part R n i(r). 



r [nm] r [nm] r [nm] 


Figure 3. Radial parts of the wave functions of the two lowest states (red and blue 
lines) for different values of Vqd: -4 meV (a), 1 meV (b), and 5 meV (c). The black 
dotted line shows the corresponding shape of the confining potential V (r) 


As already mentioned, the main advantage of the DRN is the controllability of the 
shape and the distribution of the electron wave functions. For instance, if the minimum 
of the potential of the QD part Vqd is much deeper than the potential of the QR, the 


V [meV] 












































Charge transport 


6 


electrons are located mainly in the QD and the effective size of the ground state (G) wave 
function is small. On the other hand, if the ring’s potential Vqr is much deeper electrons 
occupy mostly states in the QR part and the G wave function is much broader. What 
is more, by fine-tuning the confinement potential we can control positions of individual 
states. This way we are able to have, e.g., the ground state located in the QD, whereas 
the first excited state (E) in the QR and so on. The distributions of the wave functions 
of the two lowest energy states for three different values of Vqd are presented in figure 
3. One can see there the case where both the G and E wave functions are in the QD 
for Vqd = —4 meV (figure 3a), the E wave function is in the QR, whereas the G wave 
function is still in the QD for Vqd = 1 nieV (figure 3b), and finally, for Vqd = 5 meV 
both the G and E wave functions are in the QR (figure 3c). 

The DRN is coupled via tunnel barriers to the S and D electrodes. We assume that 
one or a few (no) single-electron states are in the bias window 

p s > 6i> hd, i = 0,..., n 0 - 1, (4) 


where fis and /i£> are the respective chemical potentials, i represents a set of quantum 
numbers (n, /) and we start the numbering of the energy levels from % — 0 (n = 0, l = 0) 
(the ground state). Charge transport through a DRN depends crucially on the coupling 
strength of its states to the electrodes what is dependent on the wave function overlap 
that enters the electron tunneling matrix element. Because we are able to control the 
shape and distribution of the DRN wave functions, we can control the overlap, which, 
in turn, allows us to control the transport properties. 

The tunneling rates T’s are calculated microscopically for each state independently. 
We follow Bardeen’s approach [48], where 

T(e) = 2 tc ^ M 2( 5 ( e - e fc) • (5) 

k 


Within the framework of this method two separate sets of states are considered: one 
solves the Schrodinger equation for the DRN and one for the electrode. It is implicitly 
assumed that the wave functions of the two subsystems are orthogonal. Enforcing this 
assumption the Bardeen tunneling matrix element in two dimensions is given by 


tk 




d%Jr) 

dx 


*nl(r) 


fopk( r ) 

dx 


5 


X=XQ 


( 6 ) 


where ^.(r) are wave functions in the electrode and H ’ n i{r) are wave functions in the 
DRN. The integral is calculated along the line x = Xq (see figure 4). The explicit form 
of tfc is given in the Appendix. 


2.1. Phonon relaxation in a DRN 

When more that one state is included in the transport process one has to take into 
account mechanisms that allow transitions between the states, i.e., relaxation to the 
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Figure 4. The geometry of the setup used to calculate the couplings T between the 
DRN and the electrodes. 


lower energy states. One of the most important and unavoidable mechanisms of 
scattering in solid state systems is the interaction with lattice phonons. For low¬ 
dimensional semiconducting systems as QDs and QRs the dominant process is the 
electron-acoustic phonon interaction. This is because the energy distance between the 
electronic states is small comparing to the optical phonon energy [33,49,50]. 

In this section we calculate phonon emission relaxation rates due to interaction 
of an electron with piezoelectric (PZ) and deformational (DF) phonons. The electron- 
phonon scattering rate due to transition from the initial state to the final state 

(r) with emission of an acoustic phonon can be calculated using Fermi’s golden rule 
(w^f = Try, 

«*-,/ = — Ei — Hu:,), (7) 

q, A 

where q is the phonon wave vector, Ef ( E ,) is the energy of the final (initial) electron 
state, hco q is the energy of a phonon and A is the polarization index. The interaction 
operator W\(q ) is given by 

W x (q) = A x (q)e-^ r , (8) 

where A\(q) is the total scattering matrix element [33] 

i^(</)i 2 = i^ptr+ £ k z («)i 2 - (9) 

A=LA,TA 

Because in this paper we focus mainly on transport properties of the DRN, we follow 
the reference [33] and use angular averaged piezoelectric coupling matrix element for 
longitudinal and transverse phonon modes. Thus, the total electron-phonon scattering 
matrix element can be written as [33,51] 

|A(9)|2 = 2yFw( fl2 ' 9 ' 2+ P- < 10 > 

where D (P) is the deformational (piezoelectric) potential constant, p is the crystal 
density, c is the sound velocity, and V is the volume. After Refs. [33,51] we assume 
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D = 2.2 x 10 18 J and P = 5.4 x 10 20 J 2 m 2 . The total relaxation rate from state n'l' 
to state nl at T = 0 can be written as 


^n'l'—ynl 


X 


1 

47 T 2 hpc 2 


4 ( D M + p ) 



(10 sin 6 



- [ ' cos 2 (V e~ izq ° cos6 dz 

dJ-d/2 v dJ 


2 


( 11 ) 


where go = (E n 'i' — E n i)/hc and R n i(r) is the radial part of the in-plane electron wave 
function ip\\(r). The integral over z is given by / (!dg 0 cos$), where 


f(x) 


71 


7 T 2 X — X 3 


sm x 


( 12 ) 


and the integral over ([>' can be expressed by the Bessel function. Finally, the relaxation 
rate can be written as 


Wn'l'-ml 


2t r 
hpc 2 


4 (d\1 + P) 



(13) 


where 


F n 'v,ni{0) = f ( ~dq 0 cos 9 


dr J\i-i>\(q 0 r sin 9)r R nl (r)R n ' V (r). 


(14) 


Apart from material constants the factors that affect this rate are the energy gap between 
states n'l' and nl, i.e., E n n' — E n i, the mutual distribution of the wave functions given by 
R n i{r) and R n 'i'{r ), and the thickness of the structure d. Function F n n> n i(0) describes 
how phonons emitted in different directions contribute to the relaxation process. The 
relaxation through phonons emitted at a given angle 9 depends on d through function / 
given in equation (12). On the one hand, since its argument is \dq$ cos 9 the relaxation 
is independent of the thickness d for phonons with wave vectors parallel to the xy plane 
(9 = 7t/2). On the other hand, there is a strong dependence of the relaxation rate on 
d for phonons emitted in the direction perpendicular to nanostructure. Figure 5 shows 
f 2 (x). One can see there that for 9 = 0 a significant contribution to the relaxation 
comes only from phonons with wavelength larger then d. 

Figure 6 shows the square of the function F given by equation (14) for d = 5 nm. 
The value of q o and the radial parts of the wave functions R n i(r) and R n 'i'{r) have 
been calculated for two different values of Vqd- In the case presented in the left panel 
(Vqd = — 6 meV) the bottom of the QD potential is much below the bottom of the 
QR potential. Therefore, the wave functions are situated mainly in the QD part of the 
DRN, where the energy level spacing is large (see figure 2). On the other hand, the 
right panel presents the case of Vqd = 6 meV where the bottom of the QD potential 
is above the bottom of the QR potential. In this situation the wave functions are 
situated mainly in the QR part of the DRN and, as it infers from figure 2, the level 
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Figure 5. Square of the function f(x) given by equation (12). 




Figure 6. Square of the function F n n^ n i(9) given by equation (14) for two different 
values of Vqd for a DRN with d = 5 nm. 


spacing is small. Figure 6 allows one to analyze the directions of phonons emission 
in two presented cases. Namely, in the left panel (Vqd = —6 meV) function F^, v nl { 6 ) 
describing the phonon emission is peaked within the range of low values of 6. This 
indicates that phonons are emitted mostly perpendicularly to the DRN, which is the 
direction of the strongest confinement [49]. In the complementary case presented in the 
right panel (Vqd = 6meV), it is not possible to point out a specific direction of emission 
(it varies with a particular transition n'V —V nl , note the difference in the scales on the 
horizontal axes). 

Figure 7 shows typical examples of how the relaxation rates are affected by changing 
the position of the bottom of the QD potential Vqd from below to above the bottom 
of the QR potential (Vqr = 0). Figure 7a shows the relaxation rate from state 
(n = 0, l = 1) to the ground state (n = 0, l = 0). For large, negative value of 
Vqd both the wave functions Too and T 0 i are positioned in the QD what gives fast 
relaxation with rates w of the order of GHz. With increasing Vqd the excited state 
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Figure 7. Illustrations of the dependence of different relaxation rates on Vqd for 
d = 5 mn. See the text for an explanation. 


moves over to the QR leading to the decrease of the overlap between the wave functions 
what, in turn, results in the decrease of w. For further increased Vqd both the ground 
and the excited state are situated in the QR, the overlap increases and the relaxation 
rate increases again. Figure 7b presents the relaxation rate from state (n — 1, l — 0) 
to the ground state. In this case for small Vqd the ground state is localized in the 
QD, while the excited state (n = 1, l — 0) is mostly positioned in the QR. Therefore, 
the overlap of the corresponding wave functions is small, what in turn results in a slow 
relaxation. With the increase of Vqd the ground state starts to move over to the QR and 
the overlap with dRo starts to rise. One observes it as a sharp increase of the relaxation 
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Vqd [meV] 

Figure 8. Relaxation rate woi->oo as a function of Vqd for different values of sample 
thickness. 


rate Wio^-oo- However, with further increase of Vqd the ground state remains in the QR, 
yet the excited state moves over to the QD and changes sign. This altogether causes the 
decrease of the relaxation rate. Figure 7c presents the relaxation rate between states 
(n — 1, 1 = 0) and (n = 0, l = 1), a transition that is a part of an indirect relaxation 
process. What is clearly visible at first glance is the abrupt decrease of w io-h-oi when the 
states cross (compare figure 2). In the range of Vqd from -2 meV to 2 meV the overlap 
between T 10 and Toi is strong. Thus, the resulting values of the relaxation rate are 
high. 

As it was already mentioned, the phonon emission at different angles is strongly 
dependent on the distribution of the wave functions (see figure 6). This, in turn affects 
the dependence of the relaxation rates on the sample thickness d. In figure 8 we present 
the relaxation rate icoi-suo as a function of Vqd- Its dependence on d is pronounced for 
low values of Vqd, when the wave functions are situated mainly in the QD region of 
the DRN. In accordance with what was reported in Ref. [33] for QD’s, we also observe 
oscillations of the relaxation rates as a function of the sample thickness d in the range of 
negative values of Vqd (yet in figure 8 we present only a part of the noted oscillations). 
In the range of low values of Vqd the level spacings are of the order of single meV’s (figure 
2). The corresponding wavelengths of emitted phonons are of the order of the sample 
thickness (couple of nm’s). This match results in the strong increase of the relaxation 
rates for such Vqd- In the complementary range of Vqd, where the bottom of the QD 
potential is above the bottom of the QR potential, the level spacings are of order of 0.1 
meV. Then, the corresponding wavelengths of emitted phonons are significantly larger 
than the sample thickness and we do not observe any strong dependence of w oi-s-oo on the 
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thickness of the structure d. Additionally, we would like to stress the strong dependence 
of the relaxation rates for all other transitions w n 'i'^. n i on d (not shown). 

3. Single electron transistor 

We are now ready to discuss the transport properties of the DRN. We evaluate the 
current given the energy spectrum and the subsequent relaxation and tunnel rates. We 
discuss below sequential tunneling current in the Coulomb blockade regime near the 
N = 0 -H- 1 transition and neglect higher order tunneling events [52], In the pure 
quantum dot case SET is switched to the conducting state when some energy level, 
shifted by gate voltage, enters the bias window. The current exhibits then the current 
peak [5]. In case of the DRN the mechanism is different. We keep one or a few states in 
the bias window and manipulate the distribution of the wave functions to get a proper 
transistor behavior. The principle of operations in is to control the energy states and 
tunnel couplings by means of gate voltages Vqd and Vq. Many investigations of transport 
behavior make use of the high tunability of the tunnel barriers by applying voltage pulses 
[52], In our approach we utilize instead the high tunability of the electron states in the 
DRN while keeping the barrier parameters constant. We assume that the bias window 
is smaller than the charging energy Ec, so that only a single electron at a time can be 
transmitted through the DRN, \ps~ fin < Ec- I 11 order to determine for what values of 
the model parameters this condition can be fulfilled we calculate Ec for the DRN. The 
interaction energy of two particles confined in the potential V (r) was calculated using 
the configuration interaction approach [53], which is an exact diagonalization method 
for solving the nonrelativistic Schrodinger equation for a multi-particle system. From 
the single-particle orbitals we constructed the basis of the Slater determinants \S tJ ). 
Then, the two-particle Hamiltonian was diagonalized and the exact eigenstates were 
found where the v-t\i eigenfunction is in the form of the linear combination of the Slater 
determinants: 

= ( 15 ) 
r 

Coefficients c ^ were calculated by the two-particle Hamiltonian diagonalization. In 
figure 2 the energy Ec of the lowest two-particle state is shown by the black dotted line 
as a function of Vqd- It can be seen there that for the analyzed range of Vqd, Eq is 
larger then the five lowest single-particle levels. 

The current through the DRN is calculated with the help of the rate equations. 
Assuming that no energy levels he in the bias window ps — p-D, the time evolution of 
the occupation probability pj of a given DRN state j can be expressed by the following 
formula: 

( no - 1 \ rao-1 / j-1 \ 

1 - ^ dipi I + ^ d i W i^jPi ~ r f + Y. djWj^j I pj. (16) 

i =0 / i—j +1 V *=0 / 

Indices i and j denote pairs of quantum numbers (n, l ) with i — 0 describing the ground 
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{ 1 if i denotes a state with / = 0, 

(17) 

2 if i denotes a state with l > 0 

is the orbital degeneracy of i-th state (E n i = E n _i). The processes that enter equation 
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Figure 9. Processes included in equation 16. The solid (dotted) S-shape line 
represents the transfer of an electron from the source electrode to state j (from state 
j to the drain electrode). The vertical solid (dotted) lines represent relaxation from 
states above state j to state j (from state j to states below state j). 


(16) are illustrated in figure 9. The first term on the r.li.s. describes the rate at which 
electrons tunnel from the S electrode. Due to the Coulomb blockade such a transfer is 
possible only when none of the DRN states is already occupied, what is ensured by the 
expression in parentheses. The second term describes relaxation to state j from higher- 
lying states. The last term describes tunneling to the D electrode and relaxations from 
state j to lower-lying states. If j = n 0 (j = 0) the second (third) sum in equation (16) 
should be omitted since there are no states from (to) which relaxation is possible. 

As we are interested in the steady-state current we put pj = 0 in equation (16). 
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Then, introducing p* = dipi, the system of equations for pj can be rewritten as 
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(18) 


In the absence of magnetic field do is always equal to 1, but we have left it in the first 
row for consistency of notation. 

The solutions can be easily found in limiting cases of a very slow relaxation and 
of a very weak couplings to the electrodes. In the former case, when we put w^j —>■ 0 
for all i and j, the solution for a symmetric coupling (Tf = Tf ) takes on a simple form 
Pi = l/(n 0 +1). In the latter case, when we put Tf, Tf — > 0 for all i, we get p* = 0 for all 
i > 0, what means that only the ground state participates in the transport. In a general 
case the system of linear equations (18) can be solved analytically, but with increasing 
no the formulas quickly become very long. With the help of the occupation probabilities 
Pi, the steady-state current can be expressed as a sum of currents If carried by electrons 
which tunnel from the S electrode to all the states that lie in the bias window. Since in 
the steady state the currents through both the barriers are equal, the total current can 
be also expressed as a sum of currents /° ut carried by electrons which tunnel from the 
DRN to the D electrode. Then, the total current can be written as 


no — 1 no—1 

1 = E = E < 1£) ) 

i=0 i=0 

where the currents carried by electrons tunneling to and from individual levels are given 
by 

/“ = er?(l - ft), /”* = erf ft. (20) 

We assume, following the experiments in Refs. [9,10], ksT « 0.01 meV. Thus, 
in our studies both the energy spacings and the bias window are much larger than the 
thermal energy and we neglect the temperature smearing out of the energy levels. Then 
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for a forward bias (/as — /id > 0) the electrons can tunnel only from the S electrode to 
the DRN and then from the DRN to the D electrode. The tunnel rates in equation (18) 
depend on the DRN geometry. We assume it so that T for the most strongly coupled 
state is equal to 5 GHz and calculate all the other couplings accordingly. That way all 
the tunnel rates are in the range of the experimentally accessible values [54], 

We start with the low bias regime where only the ground state lies in the bias 
window (n 0 = 1). In this case the total current is given by equation (20) with the 
requirement that I™ = /g ut : 


I = e 


rgrg 

rgTTg' 


( 21 ) 


We assume symmetric tunneling barriers to the S and D electrodes Tf = Tg and 
numerically calculate their values. The mechanism of the SET is as follows: if the 
bottom of the QD potential is significantly below the bottom of the QR potential, the 
ground state is localized in the QD near the center of the DRN. Since for a deep enough 
QD potential the wave function decreases almost exponentially with increasing distance 
from the QD, its overlap with the electrode’s wave functions is negligible. With an 
increase of Vqd this state moves to the outer (QR) part of the structure where it has 
much larger overlap with the states in the electrodes. This results in a strong increase 
of the current. We calculate the current as a function of the gate voltage Vqd for three 
different values of Vq. Results are presented in figure 10. One can see in this figure that 


Vn = 10 meV Vn = 20 meV Vn = 30 meV 



Figure 10. The current as a function of Vqd for three values of Vo: 10, 20 and 30 
meV. T = 500 MHz has been assumed in all the cases. 


the steepness of the /(Vqd) characteristics increases significantly with an increase of 
Vo. If Vo is small, the ground state wave function gradually moves with increasing Vqd 
towards the outer part of the DRN. This leads to a slow increase of the current. On the 
other hand, if Vq is large, there are two pronounced minima of the confining potential 
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(one in the QD and and one in the QR) and with increasing Vqd at some point the 
ground state wave function “jumps” from the QD to the QR. It results in a step-like 
/(Vqd) characteristics. This dependence is clearly seen in the insets in figure 10, where 
the differential conductance dl/dVqp is presented for different heights of the barrier Vo 
separating the DRN. 

The low bias limit, however, not always can be reached. In many cases the evolution 
of the DRN’s spectrum shows several level crossings when Vqd is changed. In such cases 
the bias ps — dD cannot be adjusted to include only the ground state for all values of 
Vqd for which the transistor effect occurs. We see in figure 2 that depending on Vqd 
two or three (or even more) states have to be accounted for. It is known that transport 
in the Coulomb blockade regime can be suppressed by the occupation of excited states 
[55]. Therefore, we analyze below how the system should be designed to get a good 
transistor behavior in the high bias regime. The dc currents are calculated as steady 
state solutions to the rate equations [equations (18-20)] for the occupation probabilities 
for all relevant states considering all the processes transferring electrons, namely the 
tunnel rates T’s and subsequent relaxation rates w' s. It turns out that the major factor 
that determines the current is the relations between the tunnel and relaxation rates. 
When the coupling to the electrodes is small compared to the corresponding relaxation 
rates, most of electrons that travel through the DRN will relax to the ground state 
before leaving the DRN by tunneling to the drain electrode. Then, the total current 
is determined mostly by the transport only through the ground state even if some of 
the excited states are in the bias window. In this case the system behaves like in the 
described above low bias limit. This situation is illustrated in figure 11. 

What we have also noted from our numerical calculations is that if the relaxation 
rates re’s are comparable to the coupling constants T’s, then only a minor influence of 
higher excited states appears. They do not change the switching characteristics but 
increase the current amplitude. 

As follows from the calculations in Sec. 2.1 the relaxation rates are of the order 
of a few GHz in the broad range of Vqd- Therefore, one can easily obtain T’s smaller 
or comparable to the corresponding relaxation rates and consequently get the desired 
transistor behavior, as shown in the inset in figure lib. 

4. Current rectifier 

Current rectification plays an important role in electron transport and is fundamental 
for the development of novel basic elements in nanoelectronics. The attempts have 
been made to build downscaled rectifiers. Coulomb blockade rectifier on a triple dot 
system [26, 56] has been recently introduced. Rectification properties of a quantum wire 
coupled asymmetrically to a quantum dot have also been studied [57]. We propose here 
a rectifier built on a DRN which utilizes different distribution of the ground and excited 
states wave functions. 

To show the idea we discuss at first a case where the two lowest states are in the 
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a) Vo = 10 meV b) Vq = 30 meV 



Figure 11. The total current along with contributions that flow through particular 
DRN’s levels for T = 500 MHz, Vo = 10 meV (a) and Vo = 30 meV (b). In both panels 
the meaning of the lines is the same. The inset in the right panel shows the total 
current in non-logarithmic scale. 


bias window and we present the analytical solutions for the currents. 



Figure 12. The scheme of a rectifier with two energy levels in the bias window for 
forward (a) and reverse (b) bias. The red dotted line represents the tunneling rate 
between the source electrode and the ground state that is assumed to be extremely 
small. The solid zig-zag line represent the relaxation w from the excited state. 


The confinement potential of the DRN is tuned so that the first excited state 
(E) has significantly larger overlap with the leads than the ground state (G). Thus, 
the same holds true for the corresponding tunneling rates, i.e, T^ <C T^. To get 
current rectification we have to break the left-right symmetry, therefore we assume 
an asymmetric coupling to the source and drain electrodes, i.e., the drain electrode is 
much closer to the DRN than the source one. This results in the following relations 
between the tunneling rates: 

r? «r? (i = G,E), (22) 

where T^. is assumed to be so small that the current between the ground state and the 
source electrode is negligible (Tg < 100 kHz [5]). One can tune the distances between the 
DRN and electrodes to ensure the above relation. Moreover, it follows from equations 
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in the Appendix that the ratio of the tunneling rates to the ground and excited states 
is independent of this distance, so that 



(23) 


what is on line with the assumptions made, e.g., in Refs. [10,54] in the interpretation 
of their experiments. 

The proposed mechanism of the rectification is the following: for forward bias 
( i ps — Pd > 0, see figure 12a) the current from the source electrode flows to the DRN 
through the excited state and then either further flows to the drain electrode through the 
excited state or the electron relaxes and leaves the DRN through the ground state. This 
is the forward direction of the rectifier. On the other hand, for reverse bias (ns~ ho < 0, 
see figure 12b) an electron can enter either the ground or the excited state. If the electron 
enters the ground state, it cannot leave the DRN because of the negligible coupling to 
the source electrode (we assume the temperature to be low enough to prevent from 
exciting the electron to the next energy level and from tunneling back to the drain 
electrode). If the electron enters the excited state, it will relax to the ground state due 
to fast relaxation (Sec. 2). Then, in either case the electron gets stuck in the ground 
state blocking the current through the excited state by means of the Coulomb blockade. 
This is the reverse direction of the rectifier. 

To describe the proposed rectifier effect in a quantitative way we calculate the 
currents using the steady state solutions of the rate equations for the occupation 
probabilities [equations (18-20)]. Taking into account only the two lowest states we 
obtain the current If for the forward bias and the reverse current Ir for the reverse 
bias: 


If = eTi 


d\ (rg + w) 


+ dlW ^1 + 


1 E 


Ir = ~ eT 


w 


+ r£ 


G 


W 


(24) 


(25) 


where w is the relaxation rate for the first excited state and d\ is its degeneracy. These 
formulas allow us to find the conditions under which DRN behaves as a rectifier, i.e., 
for which 

If I 


'R 


> 1. 


(26) 


Since Tf. <C w we get 



rf di(rg + w) 
r s 

1 G 


I'e + d\w f 1 + 

1 G 


( 27 ) 



Charge transport 


19 


One can see that the ratio of the forward and reverse currents depends crucially on 
the ratio between couplings of the ground and excited states to the S electrode and this 
is a parameter that can easily be tuned in a DRN. 

We calculate the tunneling rates P|, and Y E , Yq and the relaxation rates w 
as a function of Vqd under the assumption that the geometry of the DRN is that the 
maximal value of Y E is 5 GHz. 



Figure 13. The couplings T between an electrode and particular states as functions 
of Vqd in linear (a) and logarithmic (b) scales. Additionally, the solid red line in the 
right panel represents the resulting relaxation rate. 


Choosing Vqd = —4 meV (only two states in the bias window) for Y E /Y E = 100 we 
get from equations (24), (25) and (27) If = 60 fA, Jr = 0.9 fA, what gives \If/Ir\ = 64. 

In this reference situation we get very small currents because both the G and E 
states are placed in the inner part of the DRN (see figure 3a). To get a stronger current 
Ip we have to choose larger Vqd f° r which some excited states are placed in the outer 
part of the DRN resulting in the larger ratio of the respective tunneling rates Tf. /Tq. 
For Vqd > — 2 meV we have to consider three states in the bias window, as can be 
infered from figure 2. Below we present two examples of our numerical results: 


Vqd 

rg/rf 

If 

Ir 

\If/Ir\ 

wave functions 

1 meV 

100 

9 pA 

0.09 pA 

100 

figure 3b 

0 meV 

200 

4 pA 

0.012 pA 

320 



Table 1. Examples of model parameters and resulting currents. 

The wave functions in these two cases differ only qualitatively and therefore only 
the case of Vqd = 1 meV is presented in figure 3b. The above example illustrates that 
it is possible to design a DRN that allows one to get a substantial forward current and 
a high degree of rectification. The Coulomb blockade due to the electron stuck in the 
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ground state will prevent from transport through any of the excited states in the reverse 
direction. On the other hand, they will all participate in the transport in the forward 
direction, what increases the total current. Because Ig is not exactly equal to zero, we 
get some leakage current in the reverse direction but still substantial rectifying behavior 
occurs. 

The crucial requirement for the rectifier is the strong difference between the coupling 
to the ground state and to the excited states. This is the point where the advantage of 
the DRN over a QD is clearly visible. We performed, for comparison, calculations for 
QD with R— 70 nm, T®/T S E = 100 and have got in the most favorable case: Ip = 9 pA, 
Ir — 1-6 pA, \Ip/Ip\ = 5.7. One cannot obtain here a high degree of rectification 
because it is impossible to change the relative distribution of the ground and excited 
state wave functions in a QD. As a result, one cannot decrease the current hi without 
simultaneous decreasing of Ip. In QDs the ratio Tf/Tg. is usually below 10, whereas in 
the DRN the confining potential can be tuned to give rf,/r;| up to 10 3 . We see that the 
tunnel coupling to the reservoirs in the DRN is tunable over a much wider range than 
in QD. This, in turn, results in larger values of currents and more efficient rectification. 

5. Summary 

The fundamental requirement for future, low power consumption electronics is to 
control and manipulate single charges or spins. Quantum dots are the most popular 
and developed few-electron systems due to the relative easiness of fabrication and 
manipulation. On the other hand, the simple geometry of QDs allows modification 
of their electronic properties only to some extent, what encourages scientists to explore 
more complex systems like double or triple QDs. In this paper we performed systematic 
studies of electronic properties of a concentric dot-ring nanostructure and showed that, 
thanks to its non-trivial geometry, the structure offers unique possibilities to manipulate 
the electron wave functions. In particular, we have shown that by simple electrostatic 
gating one can move over the electron between the outer ring and the inner dot changing 
orbital relaxation by orders of magnitude and switching the character of a DRN from 
insulating to conducting. We have demonstrated that a DRN occupied by a single 
electron can be a good single electron transistor and a current rectifier with very 
high on/off ratio. The presented so called wave function engineering technique allows 
designing the properties of a system from the lowest quantum mechanical level. This is 
exactly how modern nanotechnology works and is the way to reach the limits of device 
miniaturization and to make the next step in development of quantum computers. 

In the paper we present results for a DRN build as a InGaAs structure. However, 
the idea of a nanosystem where the wave functions can be moved over to different 
spatially separated parts can be applied also to other physical systems, for example 
graphene nanostructures [58-61]. We also neglect the spin effects. Yet, exploiting them 
would allow one to use a DRN in nanospintronics. The effect of spin on the Coulomb 
blockade has been studied in single QDs [5,55,62,63] and in double QDs [29,64,65]. A 



Charge transport 


21 


spin blockade in single electron transistor in QD resulting from spin polarized leads has 
been discussed in [66]. In the case of a DRN one would be able to independly control 
spin-up and spin-down quantum states, what would lead to spin dependent tunnel rates 
T. This, in turn, would give the possibility to control a spin polarised current. 

By studying one-electron properties we have demonstrated that the structure offers 
the unique possibilities to manipulate the distribution of the wave functions what 
influences different DRN features. It also provides great flexibility in manipulating 
many-electron states, what, however, is out of the scope of this paper and will be 
presented elsewhere. 
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Appendix A. Derivation of the tunneling matrix element 

According to Bardeen’s approach [48] the tunneling matrix element is given by 




h 


2 m* 


dS- 


tp* k (r)V^ni{r) - T nZ (r)V^(r) 


what in two dimensions can be written as 

h 


tk 


2 m* 


dy 




dx 


dx 


(A.l) 


(A.2) 


X=XQ 


where T nZ (r) is the DRN’s wave function and V’fctV) is the electrode’s wave function. 
The electrodes are modelled as half-planes and their wave functions and eigenenergies 
are given by 
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The wave function of the DRN can be written as 

= Rn,(r)e u * = Rn ,(ys r +yy'"«“i. 


(A.5) 

(A.6) 
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The derivatives in equation (A.2) are given by 


d^nijr) 

dx 


dR n i(r) x „ , . ily 

Or ^ + y ^ +rUr CT* 


exp (il arctan — J 

' X' 


(A.8) 


and 


&ipk(r) 

dx 


K^-e^y sin e K ( x ~ x i)' 


(A.9) 


For each k in equation (A.l) one has to calculate k and 6 [equations (A.5) and (A.6)], 
insert equations (A.3), (A.7), (A.8), and (A.9) into (A.2) and (numerically) calculate the 
integral. It is convenient to express in polar coordinates where it takes the following 
form 
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